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Abstract. We study the averaged mean curvature flow, also called the vol- 
ume preserving mean curvature flow, in the particular setting of axisymmetric 
surfaces embedded in R 3 satisfying periodic boundary conditions. We establish 
analytic well-posedncss of the flow within the space of little-Holder continu- 
ous surfaces, given rough initial data. We also establish dynamic properties of 
equilibria, including stability, instability, and bifurcation behavior of cylinders, 
where the radius acts as a bifurcation parameter. 



1. Introduction 

The averaged mean curvature flow is a second-order geometric evolution law, 
acting on closed, compact, connected, sufficiently smooth hypersurfaces T immersed 
in R™. The evolution of T involves the (normalized) mean curvature "H(T), which is 
simply the sum of the principal curvatures on the surface. In particular, one seeks 
to find a one-parameter family of smooth, immmersed, orientable hypersurfaces 
{T(t) : t > 0} which satisfy the evolution equation 



(1.1) 



V(t)=h(T(t))-n(T(t)), t>0, 

r(o) = r , 



where V(t) denotes the normal velocity of the hypersurface T(t), Tq is a given initial 
surface, and h(T) is the integral average of the mean curvature H(T) on T. Two 



important features of the flow (1.1) are that the surface area of T(t) is decreasing 
in t and the (signed) volume enclosed by the surface T(t) is preserved, as long as 
smooth solutions exist. These features provide a starting point for our analysis of 
the dynamical properties of solutions. 



Problem (1.1) is a modification of the mean curvature flow, which is the related 
evolution equation without the integral average term h(T) in the governing equa- 
tion. These two problems have a long and rich history of investigation which we do 
not attempt to summarize here. We will simply highlight some of the important 
historical results which bear on our current investigation. 

Both the mean curvature and averaged mean curvature flows have been con- 
sidered in a wide range of settings for T. Cases considered include curves in R 2 , 
[THIHE], hypersurfaces immersed in R n , [5j Q31 [TH1 [2S] > and immersed submanifolds 
of smooth Riemannian manifolds (Ai n ,g), [T5], for arbitrary dimensions n > 2. 

The first well-posedness results for ( |1.1[ ) were established by Gage [T5] , in the set- 
ting of plane curves 7 C R 2 . Under convexity assumptions on the initial curve 70, it 
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is shown that solutions are global, convex for all time t > 0, and solutions converge 
to a circle. These results were later generalized by Huisken |18j . who proved ex- 
istence of global solutions assuming smooth uniformly convex initial hypersurfaces 
in W 1 . It is also shown in this setting that solutions converge to an n-dimensional 
sphere. Escher and Simonett |14j later proved well-posedness for rough initial data 
(in the little Holder spaces h 1+l3 ), including the possibility of nonconvex hyper- 
surfaces, and further proved that the family of n-dimensional spheres is locally 
exponentially attractive in the topology of h 1+ P . Regarding further qualitative 
properties of solutions, Mayer and Simonett [26] proved the first analytic results 
regarding hypersurfaces evolving according to (1.1) which lose embeddedness. 

Considering axisymmetric, or rotationally symmetric, surfaces, several authors 
have investigated both the mean curvature and averaged mean curvature flows in 
this setting, c.f. [5J |S] and references therein. Of particular note is a result of 
Athanassenas [5] , in the setting of smooth surfaces satisfying Neumann boundary 
conditions, regarding global existence and convergence to a cylinder for given initial 
surfaces satisfying an isoperimetric-type inequality. A recent article of Hartley |17j 
refines this result to rough initial data, also allowing for non-rotationally symmetric 
initial surfaces. An exponential convergence rate is also established, given initial 
data h 1+/3 -c\ose to a cylinder, using techniques similar to |14) . 

In the current article, we develop a function space setting adequate to establish 
well-posedness of the axisymmetric averaged mean curvature flow with rough initial 
data and to investigate geometric properties of solutions analytically. The methods 
of the paper also allow for application to higher-order nonlinear problems. For 
example the author and Simonett |23) recently applied similar methods to derive 
results for the surface diffusion flow, a fourth-order geometric evolution law, also 
in the setting of axisymmetric surfaces with periodic boundary conditions. 

We conclude the introduction with a brief outline of the article. In Section [2j we 
define appropriate spaces and introduce analytic tools that will be used throughout 
the paper. In Section [3} we take advantage of maximal regularity properties for the 
governing operator of ( |2.4[ ), in conjunction with a quasilinear structure, in order to 
establish well-posedness and regularity of solutions. We are then able to explicitly 
characterize all equilibria (stationary solutions) to the problem in Section |4j using 
results regarding closed surfaces of revolution with constant mean curvature. 

After characterizing equilibria, the remainder of the paper is dedicated to in- 
vestigating properties of two particular families of equilibria, namely the cylinders 
of radius > and families of ^-periodic unduloids. In Section [5j we develop 
nonlinear stability and instability results for cylinders, where the size of the radius 
governs the stability of the cylinder. One will see in this section how maximal 
regularity provides a crucial connection to tools in nonlinear functional analysis, 
including applications of the Implicit Function Theorem on Banach spaces. We 
conclude the article with an investigation of bifurcations of equilibria, in particular 
we use analytic methods to establish the qualitative structure of the intersections 
between the families of unduloids and the family of cylinders. 



2. Notation and Conventions 



For the remainder of the paper, we consider the case of T C R 3 an embedded 
surface which is symmetric about an axis of rotation (which we take to be the 
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x-axis, without loss of generality) and satisfies prescribed periodic boundary con- 
ditions on some fixed interval L of periodicity (we take L = [— tt, tt] and enforce 2tt 
periodicity without loss of generality). In particular, the axisymmetric surface T 
is characterized by the parametrization 



r = 



r(x) cos(0), r(x) sin(0)) : i£R, G [— t, 7r]|, 



where the function r : E — > (0, oo) is the profile function for the surface T. 
Conversely, a profile function r : K — > (0, oo) generates an axisymmetric surface 
r = r(r) via the parametrization given above. Utilizing this explicit parametriza- 
tion for axisymmetric surfaces, we can recast the averaged mean curvature flow as 
an evolution equation for the time-dependent profile functions r = r(t, x). 

Given a surface T(r), it follows that the mean curvature is H(r) = K1+K2, where 

1 , 

K\ = — -j^=^= and K2 



are the azimuthal and axial principle curvatures, respectively. Meanwhile, the 
normal velocity of T = T(t) is 

V(t) = " 



and, introducing the surface area functional 

(2.1) S(r) := f r(x)y/l + rl{x)dx, 



which measures the surface area of one period of the induced axisymmetric surface 
r(r) (modulo a factor of 2ir), the integral average of the mean curvature is 

(2.2) h(r) U{r){x)r{x)^l + rl(x)dx. 

Sir) J_„ 

Defining the operator 

(2.3) G{r):=^l + ^[h{r)-n{r) 
we arrive at the expression 



(2.4) 



r t (t, x) = G(r(t))(x), ieT,(>0 
r(0) = r 0j ieT, 



for the periodic axisymmetric averaged mean curvature flow, where T := [— 7r,7r] is 
the one-dimensional torus; with the points — ir and n identified, endowed with the 
topology generated by the metric 

dr(x, y) := min{\x - y\, 2tt - \x - y\}, x, y € T. 

There is a natural relation between functions defined on T and 27r-periodic 
functions on E, which the author explores in detail in the article [35]. Throughout 
the article, we will consider profile functions in the periodic little-Holder spaces 
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h k+a (T), which are Banach spaces for k G N := {0, 1,2,.. .}, a G (0, 1), defined by 

h k+a (T) := C^(T)"'" cfe+Q<T) and 

||r||^ := \\r\\ c * +am := £max |r«(*)| + sup 1 LI— ^ 

i.e. the closure of the smooth functions in the topology of the classic Holder func- 
tions C k+a over T. If a G R+ \ Z is fixed, then we denote by h a (T) the little- 
Holder space /iL cr J+f cr } (T), where [cj denotes the largest integer not exceeding a 
and {cr} := cr — [cj • 

In addition to functions defined on the one-dimensional torus, we consider stan- 
dard classes of regular functions between Banach spaces. In particular, given Ba- 
nach spaces E,F, and an open set U C E, we denote by C k (U,F) the class of 
fc-times (Frechet) differentiable functions mapping U into F. We also denote by 
C w (i7, F) the class of real-analytic functions, which are representable (in the topolo- 
gies of E and F) as a power series of fc-linear maps from E into F, c.f. [10 . 

2.1. Maximal Regularity. One essential tool that we use throughout the paper 
is the property of maximal regularity, also called optimal regularity in the literature. 
Maximal regularity has received considerable attention in connection with nonlinear 
parabolic partial differential equations, c.f. [IJ[2j[3l[7l[24j[28l[29]. Although maximal 
regularity can be developed in a more general setting, we will focus on the setting 
of continuous maximal regularity and direct the interested reader to the references 
[TJ [23] for a general development of the theory. 

Let /Lt G (0, 1], J := [0, T], for some T > 0, and let £ be a (real or complex) 
Banach space. Following the notation of [7] , we define spaces of continuous functions 
on J := J\ {0} with prescribed singularity at as 

BUCx-p{J,E) := ju G C{j,E) : [t i-> t^utf)] G BUC(J,E) and 

( 2 - 5 ) lim p-uMQU = o\, Me (0,1) 

||u|| Bl _„ ^supt^H^IU, 

where BUC denotes the space consisting of bounded, uniformly continuous func- 
tions. We also define the subspace 

BUCl_^(J,E) := {u G C? 1 (■/,£?) : u, u G BUC\-^{ J, E)\ , fi G (0,1) 

and we set 

BUC (J,E) := BUC{J,E), BUC^(J,E) := BUC^^E). 

If J = [0, a) for a > 0, then we set 

Ci_„(J,£) := {u G C{j,E) : u G BUCi_ /J ([0, T],E), T < sup J}, 

Ct^E) := {u G C l {lE) : u,tk G d_„(J, £)}, /i G (0, 1], 

which we equip with the natural Frechet topologies induced by Bf/C 1 _ (1 ([0,r],£) 
and BUC^^T^E), respectively. 
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If Ei and Eq are a pair of Banach spaces such that E\ is continuously embedded 
in E Q , denoted E\ ^ Eq, we set 

E (J) := BUd-^Eo), fx 6(0,1], 

ti (,! Ei(J) := BUC 1 1 _^J,E )r\BUC 1 -^J,E 1 ), 

where Ex (J) is a Banach space with the norm 

\\u\\ EliJ) ^supiWiK^lk + lKiJIk 

It follows that the trace operator 7 : Ei(J) — >• £7o, defined by 7V := u(0), is well- 
defined and we denote by 7E1 the image of 7 in Eq , which is itself a Banach space 
when equipped with the norm 

\\x\\ 7El := inf j|M| El(J) : v G Ei(J) and = x\. 

Given B G C(E\,Eq), closed as an operator on Eq, we say (Eq(J),Ei(J)J is a 
pair of maximal regularity for B and write B G MTZ^(E\, Eq), if 

j t +B,^j € A SO m(Ei(J),E (J) X7EO, MS (0,1), 

where Ci Som denotes the space of bounded linear isomorphisms. In particular, 
B G MTZ l _ l (Ei,E ) if and only if for every (/, uq) G E (J) x 7E1, there exists a 
unique solution u G Ei(J) to the inhomogeneous Cauchy problem 

(u(t) + Bu(t) = f(t), tej, 
\u(0) = Uq. 

Moreover, in the current setting, it follows that 7E1 = (Eq, ^ i.e. the trace 
space 7E1 is topologically equivalent to the noted continuous interpolation spaces 
of Da Prato and Grisvard, c.f. [TJ ED ■ 



3. Well-Posedness of (2.4) 



Well-posedness of the averaged mean curvature flow is well established in the 
literature, we reference the work of Escher and Simonett [14] for fundamental local 
well-posedness in arbitrary space dimensions with rough initial data, and we also 
reference Anathanessas 5J, who considers rotationally symmetric surfaces in arbi- 
trary space dimenions K™, satisfying Neumann boundary conditions. In the current 



periodic setting in R 3 , we establish the following well-posedness of (2.4). 

Theorem 3.1. Fix a G (0, 1) and take /1 G [1/2, 1] so that 2fi + a ^ Z. For every 
initial value ro G V M := h 2fl+a (T) D [r > 0], there exists a unique solution to (2.4), 

r(-,r ) G C 1 1 _^(J(r a ) 1 h a (T))n d_ M (J(ro), h 2+a (T)), 

on the maximal interval of existence J(tq) :— [0,i + (ro)). Further, we conclude that 
a) solutions have the additional regularity 

r(-,r )GC7"((0,i + (r ))xT) for all r eV„. 



b) (2.4) generates a real-analytic semiflow on 
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c) if there exists < M < oo so that the bounds 

• r(t,r )(x) > l/M, for all x £ T, and 

• ||r(f, ro)||/je < M, for some 6 £ (1 + a, oo) \ Z, 
hold for all t G J(r ), i/ien t + (?"o) = oo. 

To economize notation, we define the spaces, for a 6 (0, 1) fixed, 

(3.1) E :=h a (T), ^ := /i 2+Q (T), and £ M := (E ,E 1 )l tOO , /iG(0,l). 

Utilizing well-known results regarding the continuous interpolation method (•, -)° ^ 
and little-Holder spaces, we conclude 

= h 2 ^ +a {T) (up to equivalent norms), for 2ti + a </ Z, 

c.f. [22j |24]. Further, let V be the set of functions r : T -> K such that r(x) > 
for all x G T and define := VCiE^ for /x e [0, 1]. Note that is an open subset 
of Ep for all // G [0, 1]. 

Before we prove the theorem, we reformulate the problem (2.4) in order to make 
explicit the quasilinear structure of the equation. In particular, we note that 

G(r)=-A(r)r + f(r) where (A, f) : V 1/2 -> C(E 1 ,E Q ) x E , 



(3.2) 



A(r)p r- 



f(r) := 2ir 



\A + f r(x)p xx (x) 
S(r) 



Ml + r^)) 

5(r) r ' 
Lemma 3.2. Lei /Lt € [1/2,1]. T/ien 



Pxx 



(1 



pe Ei, re V 1/2 . 



(AJ)eC u [ V ll ,MK v (E 1 ,E )xE l 



G (0,1]. 



Proof. The regularity of the mappings is standard from techniques in nonlinear 
analysis, noting that every r S V M is strictly positive on T and h 1+a (T) for 

/i € [1/2, 1]. Hence, we focus on the maximal regularity result. In particular, for 
r £ fixed, it remains to show that A(r) e MTZ l/ {Ei, E ), v e (0, 1]. 

Separating the terms of the operator A(r), it follows from [2H Theorem 5.2] that 
the second-order uniformly elliptic operator with variable coefficients, 

L_* 

(l+r*) x ' 

is in the maximal regularity class M.TZ v {E-i, E ), v £ (0,1]. Meanwhile, we note 
that 



S(r) 



"(x) 



t(1 + ^)) 



Pxx(x) dx 



E 

r(x) 



<C{r) / \ Pxx (x)\dx<2'KC{r)\\p\\c^ ) 



Sir) h(l + rl{x)) 



Pxx(x) dx 



e£(h 2 +( 3 (T),E ), 



for any < f3 < a, r £ V\n-, where C(r) > depends on the /i 1+Q (T)-norm of 
r and min^gx r{x). Thus, from the interpolation inequality, in conjunction with 



the identity h 2 +^(T) = [E^E^^ 9 = 1 



a-P 



, and an application of Young's 



inequality, we realize this term as a lower order perturbation, c.f. [14, p. 2793]. 
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The operator A(r) then satisfies maximal regularity by a well-known perturbation 
result, c.f. Ullll]. □ 



Proof of Theorem \3.1\ For (i E [1/2,1), the main statement of the theorem, and 
parts (b) and (c), follow from the results [7j of Clement and Simonett regarding well- 
posedness of quasilinear problems, in the presence of sufficiently regular mappings A 



and /, and maximal regularity properties, as established in Lemma 3.2 Meanwhile, 



when /i = 1, we must consider (2.4) as a fully nonlinear problem, for which DaPrato 
and Grisvard [5] have established well-posedness and Angenent [5] has established 
scmiflow properties in the presence of maximal regularity properties for the Frechet 
derivatives DG(p). Indeed, computing the Frechet derivative DG(p), one sees the 
structure 



DG[r)p = ~ ^wr L ori^ dx + B{r)p ' 

for p £ Ex, where B(r) is a first-order differential operator. Hence, we apply the 



same perturbation technique as in the proof of Lemma 3.2 to establish maximal 
regularity properties for DG(r), r £ V\. 

To prove the additional regularity of solutions, part (a) of the theorem, we rely 
on a paramater-trick that goes back to Masuda [25 and Angenent [31 0] , where 
one can introduce parameters and apply the implicit function theorem to obtain 
regularity results for solutions, see also O H3J [23] • In particular, we define the 
translation operator T a : T — > T which takes x E T to the unique element of T 
which resides in the coset [x + a] E R/27rZ, for adl. Given a solution r(-,ro) to 



(2.4), for some r E Vu, we take t\ E (0,t + (r o )) and consider the function 

r Ka (t, x) := r((l + A)i, r )(T ta x), t E I := [0, tj, x E T, 

for (A, a) E (—5, 5) 2 , with S > chosen sufficiently small. It follows that r\ :a E Ei(I) 
and is a zero of the operator $ : (E x (7) n C(I, V)) x (-5, 5) 2 -> E (7) x E^, 

$(v, (A, a)) = (d t v - (1 + A)G(u) - ad x v, - r ), 



where the spaces Ej(J) are defined as in (2.6) above. One quickly verifies that $ is 



real-analytic on the domain specified and further, the Frechet derivative of $ at r 
satisfies 

A$(r,(0,0)) = (d t -DG(r), 1 )e£ isom (MI)>Ml) X ^) 

by way of maximal regularity, c.f. [71 Theorem 6.1]. The result now follows by an 
application of the implicit function theorem, c.f. [2"5l Theorem 2.4]. □ 



4. General Properties of (2.4) 

We move on to investigate geometric properties of solutions. To begin, we verify 
that two important features of the averaged mean curvature flow hold for solutions 
of (2.4), namely volume-preservation and surface-area reduction. Following these 
verifications, we characterize the collection of equilibria to (2.4), which will be the 
central focus of the remainder of the article. 
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4.1. Volume Preservation and Surface Area Reduction. Suppose that r <E 
C\_ (J(ro),E ) n Ci_ M (J(ro), E{) is a solution to (2.4) for some initial value r g 



[1/2, 1]. Looking at the evolution of the surface-area functional S, making 
judiscious use of integration by parts and periodicity to cancel boundary terms, we 
see that 

|s(r) = / (^W) + ^h^d x ) G(r)dx 
at Jt \ \Jl + r x (x) J 

H(r)(x) \h(r)-H(r)(x)]r(x)y/l + r2(x) dx 



h(r) — %{r){x) r(x)\J\ + r x (x) dx, 



where the last equation follows from 

/ h(r) \h(r) - H(r)(x)]r(x)y/l + r%(x) dx 
Jt L ^ 



h(r) / r(x)y/l + r%(x) dx) - h(r) / H{r)(x)r(x)^/l + r 2 x (x) dx 



= S(r) h(r) - S(r) h{r) 



= 0. 



Thus, it follows that the surface area of the induced surfaces of revolution T(r(t)) 
is decreasing in t. In fact, we can see that the surface area is constant if and only 
if the mean curvature is constant, which coincides with an equilibrium surface. 
Meanwhile, looking at the evolution of the volume functional 



F(r) 



r 2 (x) dx, 



we see that 
d_ 
dt 



F(r) = 2 1 r{x)G{r){x) dx = 2 / r{x)^l + r%{x){h{r) - H{r){x)) dx = 
Jt Jt 

Therefore, taking into account the regularity of the operator F in the topology of 
E^, it follows that F{r{t)) = F(r ) for all t G J(r ). 

4.2. Characterization of Equilibria. Considering the equilibria, i.e. the steady 
states, of the problem (2.4), it follows immediately that f € is an equilibrium 



if and only if G(f) = 0. Moreover, from the structure of the equation we conclude 
G(r) = if and only if T-L{f) = h(f), i.e. the mean curvature of the surface 
of revolution r(f) must be constant on T. Using results of Delaunay [TT] and 
Kenmotsu [20] regarding surfaces of revolution with constant mean curvature, we 
characterize all equilibria of (|2.4|) as follows. 



Theorem 4.1 (Characterization of Equilibria). Fix a e (0, 1). Then f e h 2+a (T) 
is an equilibrium of the flow ( |2.4[ ) if and only if f is a 2n-periodic undulary curve. 
Moreover, up to a translation along the axis of rotation, r can be expressed explicitly 
by the parametric equation (with respect to the arclength parameter s), 



R(s;H,B) 



l + Bsm(Ht) 



L/2H V 1 + B 2 + 2B sinCHi) 



dt 



^/\ + B 2 + 2Bsm(Hs) 
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for constants H > and B £ (—1,1) which satisfy the relationship 
tA^ nn f 3 "^ 1 + Bsint 

(4.1) = / dt , for some keN. 

k J n/2 Vl + B 2 + 2Bsint 

Proof. Delaunay proved that every closed surface of revolution in M 3 with 
constant mean curvature is either a sphere, a catenoid, a nodoid or an unduloid. 
Looking at the associated profile curves, the only surfaces which fit into our current 
setting (taking into account periodicity, embeddedness and regularity constraints) 
are the family of 27r-periodic unduloids, for which Kenmotsu [20] derived an explicit 
parametrization of profile curves. The parametrization presented herein is a slightly 
modified version of Kenmotsu's formula, where we have adjusted the graph to be 
symmetric about the j/-axis and to exhibit ^-periodicity in the x-variable, for a 
given value k € N. □ 

To stimulate the reader's curiosity, we include here some graphs of undulary 



curves r which, by the previous theorem, are the only equilibria of (2.4) (up to 
shifts along the x-axis). 




Figure 1. Families of 2ir periodic undulary curves (with k = 1) 
with values between B = — .99 and B = 0.99, as indicated. 



5. Stability Behavior of Cylinders 

For the remainder of the paper, we equate the constant function r*(af) = r* > 
with the cylinder r(rv) with radius r* and length 2-7T. The stability of cylinders 
under the averaged mean curvature flow is a particularly vibrant question, in ro- 
tationally symmetric settings. Considering n-dimensional smooth surfaces of rev- 
olution with Neumann boundary conditions, Athanassenas [5] showed that initial 
surfaces satisfying a particular perimetric-type inequality, comparing the surface 
area, enclosed volume, and length of the interval of revolution, will have global solu- 
tions which eventually converge to an n-dimensional cylinder. More recent results 
by Hartley [T7] indicate that one may be able to generalize this convergence result 
to rough initial data without requiring rotational symmetry of initial data, though 
we again see a perimetric-type inequality in his results. Athanassenas notes in [5] 
that the volume constraint imposed by the perimetric-type inequality eliminates 
the possibility of unduloid equilibria, in the absence of which one can concludes 
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convergence to cylinders. However, in the current paper, we are interested pre- 
cisely in the family of unduloids and how they interact (bifurcate from) the family 
of cylinders, for which a finer analysis is necessary. 

In this section, we will demonstrate that the family of cylinders with radius 



r* > 1 is exponentially stable, c.f. Theorem 5.3 Moreover, we establish instability 
of cylinders with radius < r* < 1 and develop a setting within which dynamical 
systems techniques become accessible. In the next section, we will show how this 
setting can be applied to study the bifurcation which occurs at the critical radius 
r* — 1, and subsequent bifurcations at smaller radii, c.f. Theorem |6.1| 

We first note that functions in the little-Holder spaces, r 6 ^(T), obey the 
Fourier series representation (c.f. [22] ) 

r(x) = y^f(k) e ikx , i€T; where f(k) := [ r(x)e~ ikx dx. 

* — ' 2-7T If 

We then compute the linearization 

(5.1) DG(r*)p=-DH(r it )p+~ ( [DU(r+)p](x) dx, P eh 2+a {T) 
(5-2) = (p r~ 2 + p xx ) - — f_{p(x) r~ 2 + p xx {x)) dx, 



I 



and realize it in the form 

[DG(r*)p](x) = Yl (r^ 2 -k 2 )p(k)e lhx , peh 2+a (T), xeT. 
fcez\{o} 

From this expression, it is easy to verify that the point spectrum of DG(r t ) - which 
must in fact coincide with the entire spectrum a(DG(r+)) due to compactness of 
the embedding h 2+a (T) ^ h a (T) - is 

(5.3) a(DG(r+)) = {0} U {(r~ 2 -k 2 ):keZ\ {0}} . 

Notice that, for r* > 1, this spectrum is contained in the left-half of the complex 
plane, which we can take advantage of in order to establish stability of the equilib- 
rium r*. However, the presence of in the spectrum is troublesome, as it indicates 
the presence of a center manifold which disturbs, but does not derail, our stability 
argument. To get around this hurdle, we proceed to reinterpret the operator G, and 
subsequently DG, in a reduced setting where this eigenvalue no longer shows up. 



The following technique takes advantage of the volume-preserving nature of (2.4) 
in order to investigate local properties of the problem, it is motivated by methods 
of Prokert [37] and Vondenhoff [3D] and employed by the author in [33] . 

5.1. Zero Mean Functions and the Reduced Problem. For this subsection, 
we consider an arbitrary > and a £ M. + \ Z, unless otherwise stated. We first 
introduce the mappings 

P r := r - — [ r(x) dx, and Q := 1 - P 

which define projections on h a (T). Denote by h^T) the image Po(h a (T)) , which 
coincides with the zero- mean functions on T in the regularity class h a (T) . We then 
have the topological decomposition 

h a {T) = /^(T)®Q (^(T)) =/i£(T)®M. 
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In the sequel, we equate the constant function [t](x) = rj\ 6 Q^ih 17 (T)) with the 
value 77 € R, and we denote each simply as 77. 
Consider the operator 

(r,f, 77) ^PQi^ — r, F(r) — F(r))J, with F(r) :— / r 2 (x) dx, 



<I> 

which is a real-analytic operator from h a (T) x h$(T) x K into /iq(T) x R. Notice 
that $(rv,0,r*) = (0,0) and 

D x $(r*,0,r*) = (p , 2 / r+id h <, (T) dx 

(5.4) v 

= (P 0) 47rr,Q ) G £. lsom (fc CT (T), ^(T) x M), 

i.e. the Frechet derivative of $ with respect to the first variable, at (r*, 0, r*), is a 
bounded linear isomorphism. Hence, it follows from the implicit function theorem, 
c.f. [TU1 Theorem 15.3], that there exist neighborhoods 

(5.5) (0, r,)eU = U xU 1 C h$(T) xR, e U 2 C h a {T), 
and a C u function ip* : U — >• U2 such that, for all (r, f, 77) G J7 2 x [7, 

$(r, f, 77) = (0, 0) if and only if r = ip*(r,r)). 

Note that the quantity F(r) corresponds to the volume over T enclosed by the 
surface T(r) (modulo a factor of 2n), which is a preserved quantity for (2.4); i.e. 
if r(-,ro) is a solution to (2.4), for an appropriate initial value r$ G /i 2m+q (T), it 
follows that F(r(t)) — F(r ) for all t G J(ro). Then, we can interpret t/>*(-,?7) as 
a lifting operator which takes a zero-mean function f G hg(T) to an associated 
function in the equivolume set 

M; :={reh°(T):F(r)=F(r,)}, 

which contains all the profile functions corresponding to surfaces enclosing the same 
volume as the cylinder 77. In particular, if we choose 77 so that F(ro) — F(r}), it 
follows that r(t, r ) G M™ for all t G J(r ). 

Remarks 5.1. For > fixed, we can immediately state the following properties 
of ip* which come from the definition and elucidate the relationship between Pq and 
A- 

a) Potp±(r,i]) — f for all (f, rj) G [/. 

b) Given r G ?/>*([/) n M°, it follows that ip*(P r, 77) = r . 

c) i^*{0, 77) = 77, for 77 G Ux. This and the preceding remark follow from the fact 
that F(rj) is injective as a function on (0, 00). 

d) It follows from the identity $(^*(f, 77), f, 77) = (0,0) and differentiating with 
respect to f that 

Di^*(0,Jj)fc = fc, /i G ftJ(T), 77 G C/x. 

e) Note i/j*{Uo,ti) c .M^ for 77 G t/x and we have the representation 

Hence, we see that D E/2 is the graph of a M-valued function over Hq (T) 
which is therefore a Banach manifold, with real-analytic local parametriza- 
tion ^0,77). 
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f) A priori, 77) depends upon the parameter a, however it follows from the 
preceding representation that 

M;v) ■ u nt4(T) -> h*(T), aeR+\Z, 

so that ip* preserves the spacial regularity of functions regardless of the 
regularity parameter a with which ip* was constructed. However, notice 
that the neighborhood Uq will remain intrinsically linked with the parameter 
which was used to construct ip*. 

Fix r* > and define the reduced governing operator 

(5.6) g(f, V ) := P G(Mr,v)) 

which acts on (f , 77) in a neighborhood of U, as constructed for ip* above. We need 
to be careful with the regularity of the zero-mean functions f that we plug into 
Q, so we will assume throughout that a € (0, 1) is given, and then we define the 
lifting ip* and the related neighborhood U — Uq x U\ within the setting of /iq(T). 
Then, we take Q{-,rj) as acting on the functions f 6 Uq h^ +a (T), r\ € U\. With 
this reduced operator, we defined the ^-dependent reduced problem 



(5.7) 



f t (t,x)=g(f(t,x),7]), f>0,ieT, 
r(0,a;) = (Poro)(aO, x e T. 



Theoretically, we can choose the parameter rj g U% arbitrarily, however, in practice, 
we will choose 77 for which F{rf) = F(tq), which essentially gives us the freedom to 
consider non-volume preserving perturbations ro of the cylinder r* in our stability 
analysis. 

Fix a € (0, 1) and we denote the spaces 

F :=h%(T), F 1 :=^ +C (T) > and F„ := (F^F^, p 6(0,1), 

so that F^ = PqE^ for /i e [0,1]. Moreover, for /i G (0,1] and closed intervals 
J C R + with 6 J, define the spaces 

E (J) :=BUCi-n(J,Eo), 

E x ( J) := BUCl_^(J, E ) n BUd-^J, Ej), 

and 

F (J) :=BUCi_„(J,F ), 

Fr(J) := BUCl_^J,F ) n BUd-^Fi), 



within which we will discuss solutions to (2.4) and the reduced problem (5.7), 
respectively. We can immediately see how the lifting operator t/i* connects these 
two problems. 

Lemma 5.2. Fia; r t >0,r/e[/i and J := [0, T] /or T > 0. TTien 

(5.8) M-,V) :¥ 1 (J)nC(J,U ) — >E X (J), toift ^(f , ^(t) := ^(f (i), 77). 

Moreover, if ro := -Po r o €E Pu and f = f(-,fo) € Fi(J) D C(J,Uq) is a solution to 
(5.7), /or some /1 € [1/2, 1], then r :— ^*(f, rj) is the unique solution on the interval 
J to (2.4), with initial data r :— ip+ifo, 77) G 
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Proof. Because of the local nature of the operator we can assume, without loss 
of generality, that the neighborhood U = Uq X U\ is chosen sufficiently small to 
ensure that -0* is in the regularity class C u and the bounds 

(5-9) \\i>(r,v)\\E <N and ||£>iV>(f, v)\\c(f ,e ) < N 

hold uniformly for (f, rf) £ U. Further, the embeddings 

(5.10) Fi(J)-> BUC(J,F^^ BUC(J,F ), ^ [1/2,1], 

follow from [lj Theorem III. 2. 3. 3] and the continuous embedding of littlc-Hbldcr 
spaces, respectively. It is then straightforward, utilizing the representation for tp+ 
given in Remarks 5.1 ^e) and these bounds, to verify (5.8). 

Meanwhile, by (5.10) we have ro := ipi<(fo,ri) £ V^, and so it follows from 
Theorem 3.1 that there exists a unique maximal solution 

r(;r ) G C\_ J J(r ), E ) n Ci_ M (J(r ), 



to (2.4). Define p(-) := ip(r(-),r]) as indicated and it suffices to show that pt(t) = 
G{p{t)) for t € J := (0, T], since this will imply that p{t) — r(i,r ) by uniqueness 
and maximality of the solution r(-,fo). So, let t € J and consider the auxiliary 
problem 

f 7 (r) = G( 7 (r)), forre[0, £ ], 
\l(0)=p(t), 



which has a unique solution 7 e C 1 ([0, e], i?o) H C([0, £•], by Theorem 3.1 



pro- 
vided we choose £ > sufficiently small for the particular value p(t) £ E\. By 
Remarks 1 5 . 1 1 and preservation of volume for solutions of (2.4), we have the represen- 
tation j(t) — tp i ,{P 'y{T),vi) 1 t 6 [0,e]. We then conclude the proof by computing 

G(p(t)) = 7 (0) = d T (MPol(r), V)) = C 1 ^(Po7(0), v)Poj(0) 

T = 

(5.11) = D 1 ^(P p(t),r 1 )P G(p(t)) = Di^(f (t), rj)g(f{t),n) 

Considering the stability of r*, we compute the linearization 

DiGif, V ) = P DG(Mr, v))DiMr, »?), (f , f) G tf, 
which simplifies at r — 0, using ( |5.1| ) and Remarks |5.1[ d), 

(5.12) D x g{Q,rf)p = P (DG(r))^D 1 ^(0,ri)p = -P (p DH(v))p = -P a DH( V )p 



□ 



for p £ /iQ +a (T). Recalling (5.3), and the reduced domain of definition for £/, we 
conclude 



(5.13) 



aiD^O, 77)) = {(?T 2 ~ k 2 ) : k £ 7L \ {0}} 



5.2. Stability of Cylinders with > 1. Define the exponentially weighted 
maximal regularity spaces 

F,-(K + ,a;):={/:(0,cx))^Fo | [t ^ e ut f(t)] £ F J -(K+)} ) weR,j = 0,l, 

which are Banach spaces when equipped with norms \\u\\f. (r + , w ) := ||e LJt u||F :J (R + ). 
We use these spaces in order to show the following stability result, which essentially 
tells us that /i 2/i+a -small perturbations of the cylinder > 1 will have global 
solutions which converge exponentially fast to a cylinder rj, which is close to r*. 
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Theorem 5.3. Fix a G (0, 1), /i G [1/2, 1] so that 2/i + a ^ Z, and /ei r* > 1. 

There exist positive constants e = £(?'*), <5 = <S(f*) did a; = w(r*,<5), smc/i i/ia£ 
problem (2.4 1 wit/i initial data ro G B^jt+a (r*, e) /ias a unique global solution 

r(;r ) G C 1 1 _ M (K+, ft Q (T)) n Ci_ M (R+, n 2+Q (T)), 

and i/iere exist 7/ = 77(7*0) G (r* — 5, r* + 5) and M = M(a) > /or which the bound 

^-"Hrftro) - r?|U=+« + ||r(t,r ) - »?|U"m+- < M e -""||r - r*|| h 2„+ Q 

ZioZds uniformly for t > 0. 

Proof. We demonstrate this result by first showing that is exponentially stable 
in the reduced problem (5.7 1, and then we lift solutions using ip* and show that 
exponential convergence is preserved in the lifting process. 

(i) Fix S G (0, r* — 1). Notice that the linearization of G(0,il) nas the structure 

DiQ(0,r))p = Pa(h xx + -rf 2 pj = (h xx + rf 2 - J (p xx (x) +r)~ 2 p(xfj dx, 

from which, similar to the argument presented in Theorem 3.1 we realize DiG(0, rf) 
as a lower order perturbation of the second order derivative operator d 2 . Hence, 
it follows that DiG(0,r/) G H(Ei,Eq); the class of infinitesimal generators of ana- 
lytic semigroups on Eq with domain E\. Moreover, it follows from (5.13) that the 
spectral type 

typeiDtQiO, rj)) < 1 T ^* ~JJ < 0, for all n G (r* — 5,r± + 5). 

Hence, if we choose oj G (o, ^(r'-g)* 1 ) , then ^ follows by [l :i Theore m III. 3. 4.1 and 
Remarks 3.4.2(b)] that the exponentially weighted spaces , 

Fo(R+,w)> FiCK+.w)), 

are a pair of maximal regularity for D\Q(0, 77), for all 77 G (r* — 6, + 5). 
Define the operator 

K(r,r ,r)) := (d t f -G(r,r)), r yf-f \, 

acting on U := ^Fi(M+, w)nC(M+, U )j x (UqCiF^ x U\ , which is an open subset 
of Fi(E+, u) xF,,xl. It follows that 



/CGC"(u,F (M + , W )xF M ), 



by analyzing the mapping properties of the individual operators 7, dt and Q on 
their associated domains of definition. The fact that G(-,T]) maps (Fi(R+,w) n 



C(R+, Uq)J into Fo(M+, w) follows by utilizing the representation (3.2 ), and bound- 
ing individual terms of the resulting expression using the embeddings (5.101 and 
the exponential boundedness of functions in Fi(R + , w). 
Meanwhile, notice that /C(0, 0, 77) = (0,0) and 

^1^(0,0,77) = (d t -D 1 g(Q,ri),^ G £ isom (F 1 (R+,u;),Fo(R + ,w) x F^, 
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for r\ € (r* — 6, r* + 5). Thus, by the implicit function theorem there exists an open 
neighborhood (0,r*) G U C P M x K and a C" mapping k : U -> Fi(R+,u>) such 
that 

/C(«(fo,J7),ro,r?) = (0,0) for all (r ,r})eU. 



In particular, K(fo,rj) is a global solution to (5.7 1 with parameter 77 and initial data 
fo 6 which converges exponentially fast to 0. Without loss of generality, we 
assume U C U\ x (r* — <5, 7'* + <5). 

(ii) Choose e > so that, for every r £ BB jl (r + ,e), there exists (r, 77) G J7 
for which $(r ,r, 77) = (0,0). Let r g Bjs (r*,e), fix 77 G (7%. — 5, + 5) so that 
P( r o) — an d define 



Using Lemma 
the unique glo 



5.2 



one verifies that r G Cf_ 



_,/Y*(T))nCi_ M (R + ,/i 2+Q (T)) is 



Dal solution to (2.4) with initial data j-q. Hence, it remains to show 



exponential convergence of r to the cylinder 77. Noting that k(0, 77) = 0, and utilizing 
Remarks 5.1 and the mean value theorem, we compute 

r(t)- v = (P Q + Q ) (^(s(P r , l)(t), v) - M*(P, v)(t),v) 

= n(P r , rj)(t) + Q f^*(«(fl)r 0) »?)(*), 77) - ^*(/s(0, »?)(*),»?) 



«(Por , rj)(t) 



1 

2^ 



Pit/v (TK(P r , 77) (t) , 77) K(P r , 77) (rj, x) dr efcc, 



T JO 



for t > 0. We bound the terms 

e w V-"|KP ro,r7)(t)|| Fo 



and 



I K (P r , 77) (t) || 



uniformly for t > using the property «;(Poro,7y) G Fi(K + ,w) and (5.10), respec- 
tively. Meanwhile, bounding integral terms in the C(T)-topology and using the 
bounds (5.9), we get 

(5.14) e^V-lK*) - T7|k < (l + ciJv)||«(Poro,r7)||F 1 (R + ,a,), * > 0, 



(5.15) e wt ||r(t)-77|| B(i < ( C2 + C3 7v)||«(Por ,77)|| Fl(K+itj) , t > 0. 

Here the constant ci comes from the embedding Pi <^-> Po , and the constants C2 and 
C3 come from the embeddings (5.10). Finally, by the regularity of k, we may assume 
that U was chosen sufficiently small to ensure that D\n is uniformly bounded from 
U into Fi(R + ,w). Recalling that k(0, 77) = 0, we have 



(5.16) 



||«;(Poro,r/)|| Fl(R+jW ) < / |pi«(rPoro,»7)Poro|| Fl(K+iU)) dr 

<J 

< M\\P r \\ Fu < M\\r -rJ\ Eu , 



where M := \\P \\ sup^ ^ \\Di n (r, 77)| U(f )> ,f 1 ( 
(5.16) and the inequalities (5.14)— (5.15). 



-,")) 



. The claim now follows from 

□ 
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5.3. Instability of Cylinders with < r* < 1. 

Theorem 5.4. Let r* 6 (0, 1) and p 6 [1/2, 1] be fixed with 2p + a ^ Z. TTien £/ie 
equilibrium r+ of (2.4 1 is unstable in the topology o//i 2 ' I+Q (T) /or initial values in 
h 2 ^+ a (T). 

Proof, (i) Let r* 6 (0,1) be fixed, and let L := DG(r+). We can restate the 
evolution equation \2A\ in the following equivalent form 

( 5 17 j ipt~Lp = g(p), t>0 

where g(p) := G(p + r ir ) — hp. Using the quasilinear structure of G it is not difficult 
to see that for every j3 > there exists a number eo = £o(/3) > so that 

(5.18) \\g(p)\\ Eo </3||p|k, forall p eB B „(0,e ) n£i. 

It follows from ( |5.3| ) that 

o-(i) n [Rcz > 0] ^ 0, 

and we may choose numbers lj, 7 > such that 

[w - 7 < Rez < u> + 7] n o-(X) = and rr + := [Re z > lo + 7] n <r(L) ^ . 

Define P + to be the spectral projection, in Eq, with respect to the spectral set 
er + , and let P_ := 1 — P + . Then P + (Eq) is finite dimensional and the topological 
decomposition 

E = P+(E Q ) (B P-(E Q ) 

reduces L. Hence, L = L + 0L_, where L± is the part of L in P±(Eo), respectively, 
with the domains D(L±) = P±(Ei). Moreover, P± decomposes E\ and, without 
loss of generality, we can take the norm on E\ so that \\v\\e 1 = \\P+v\\e 1 + \\P-v\\e 1 - 
Note that 

cr(i_) C [Rez < oj - 7], v(L+) = cr + C [Rez > w + 7], 

which implies there is a constant Mo > 1 such that 

l|e-P-IU (a) <M 0e <-* 

lle-'+'P+IUia) <M„e-'"+», t>0 

where {e L ~* : t > 0} is the analytic semigroup in P_(Eq) generated by L_ and 
{e i+t : i G IR} is the group in P + (Eo) generated by the bounded operator L + . 

From [221 Theorem 5.2] one sees that (E (J),Ei(J)) is a pair of maximal reg- 
ularity for —L and it is easy to see that — L_ inherits the property of maxi- 
mal regularity. In particular, the pair (P_(E (J)), P_(Ei(,J))) is a pair of max- 
imal regularity for — L_. In fact, since type(—uj + L_) < —7 < we see that 
(P_(Eo(]R+)), P_(Ei(IR + ))) is a pair of maximal regularity for (lj — L-). This, in 
turn, implies the a priori estimate 

(5.20) Ue-^HkcM < M i(lKlk + l|e-"7l| Eo (J r )) 

for Jt := [0, T], any T 6 (0, 00) (or Jt — M+ for T — 00), with a universal constant 
Mi > 0, where w is a solution of the linear Cauchy problem 

( W ~ L ~ W = L with (j> ) G (c((0,T),P_£ o ),P-^o). 

u>(0) = w , V / 
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(ii) By way of contradiction, suppose that is stable for (2.4 1. Then for every 
e > there exists a number S > such that (5.17) admits for each r £ Me (r+, 8) 
a global solution r — r(-,r ) £ Cl_^(R + , E ) n Ci_ M (M+, which satisfies 

(5.21) IKt)K<e, *>0. 

We can assume without loss of generality that j3 and e are chosen such that 

(5.22) 2Co(M + A/i 7 )/3 < 7 and e<e (J3), 

where Co '■= Tnax{\\ P^\\c(e ), \\P+\\c(e )}- As P + (Eq) is finite dimensional, we may 
also assume that 

\\P+v\\ E „ = \\P+v\\ Eo , v£E , v€{n,l}, 
using the fact that P+E C D(L n ) for every n £ N, c.f. [231 Proposition A. 1.2]. 
CLAIM 1: For any initial value r$ £ B.g (r*, 5), P+r admits the representation 

/oo 
e L + {t - s) P + g{r{s))ds t>0. 

For this we first establish that, for ro £ 1^^, 5), 
e- u VeBCi_,(l + ,£ 1 ) : 



u £ C((0,oo),-Ei) : sup ^""IK*)!!^ < oo 

tGR+ 



The mapping property g : Ei(Jy) — > E (Jr) follows analogously to the mapping 
property derived for Q in the proof of Theorem |5.3| above, < T < oo. Together 
with the inequalities (5.18) and (5.20), this yields 



(5.24) 



\\e~^P-r\\ Bl _ AJ ^ El) < M 1 (||P_(r - r*)lk 

+ C p\\e-^P + r\\ Bl _ MTtEl) + C /3||e— *P_r|| Bl _ (l(jT)Bl) ) 



for any < T < oo. Due to (5.22), we have MiCqP < 1/2 and can further conclude 
(5.25) 



e-^P-rWs^j^ < 2M 1 (\\P_(r -r*)\\E tl +C p\\e-" t P + r\\ Bl _ AJTtEl) 



It follows from (5.21) that 

t i-P|| e -<- *P +r (t)\\ El < ^-^e-^ColkC*)!!^ <C C l£ 

where C\ :— sup{t 1 ~' 1 e~ wt : t > 0} < oo. Inserting this result into ( 5.25[ ) yields 

(5.26) \\e- ut r\\ Bl _^ jT>El) < 2M 1 ||P_(r - r*)lk + (2M 1 C /3 + l)C Q C x e < C 2 

for any < T < oo. However, since T is arbitrary and ( 5.26[ ) is independent of 
T we conclude that e~ wt r £ PCi_ M (R+, E x ), for any initial value ro £ M B (r+,S). 
Next we note that, for s > t, by (5.19) 



(5.27) 



|| e £ + (t- S ) P+g(r(s)) | Un < Af C /3 e (w+ ^ (t - s) ||r( S )|| £l 

<MoC7 /9e w *e^*-)^- 1 ||e-'-r||fl l _ (l(Rf , BO , 



by which the integral in (5.23) exists for t > 0, convergence in £i. Moreover, 

<e^M C CMe-^r\\ Bl _^^ El) , 

Eq 



(5.28) 



/OO 
e i+(*-*)p +5 ( r ( a )) ds 



IS 
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where C 3 := sup { / t °° e 7 '* s 's M 1 ds : t > 0} < oo. Noting that w — P + r solves 
the Cauchy problem 

\w - L + w = P+g{r), 
\w(Q)=P + (r -r*), 
it follows from the variation of parameters formula that, for t > and r > 0, 



P+r(t) = e 



P+r(r)+ / e L + {t - s) P + g{r{s))ds. 



This representation holds for any r > and the claim follows from ( |5 . 19 1 and (5.21 ) 
by sending r to oo. 

CLAIM 2: If r E M E)i (r*,6) and ||r(t)|| Bf< < e for all f > 0, then 

||P+(r - rv)|| E „ < 2Af Af 1 C7 3 ||P_(r - r,)||^. 



From ( |5.23[ ) and |5_27f follows 

\\e- ut P + r 



(5.29) 



Si _ M (R + ,E ) 



< 



7 



|e-"*P + r|| B . 



le-^P-rllB,,^ 



,Ei. 



where we have used the fact that sup t>0 {t 1 M J t °° e 7 ^' ^s' 1 1 ds} < 1/7. Adding 
the estimates in (5.24 1 and (5.29) and employing ( |5.22 1 yields 

(5.30) 



l|e" 



,E X ) <2M 1 ||P_(r -r*)|U (j 



The representation (5.231 in conjunction with (5.28) and (5.30) then implies 
(5.31) 

||P+(ro ~ tv)||e„ < M C C 3 /3\\e- ut r\\ Bl _^ +}El) < M C 3 \\P.(r - r*)|U M , 

where the last inequality follows from 2CoMif3 < 1. We have thus demonstrated the 
claim, and the theorem follows by way of contradiction. In particular, note that if 
tq E M E ^ (r*, 8) is chosen with ||P_(ro — r*)|| E — 0, then it must hold that ||P + (ro — 
r *)ll-EV, — 0, so that r = r*, which contradicts the stability assumption. □ 

6. Bifurcation From Cylinders 
We conclude by investigating the interactions between the family of cylinders, 



which we will informally consider as the trivial equilibria of (2.4), and the families 



of undulary curves, the non-trivial equilibria. In particular, by restricting (5.7) 



to a problem on profile functions which are symmetric about the y-axis (i.e. even 
functions on T), we observe subcritical bifurcations from the family of cylinders 
which occur at the cylinder of radius for any £ E N, where the reciprocal of 
the radius r* is taken as a bifurcation parameter. Most of the terminology and 
notations employed in this section coincide with those of Kielhofer [3T] . 
To begin, define the operator 

G(r,X) := P G(Mr,r*)), A := 1/r*, r* > 0, 

and we immediately note that 

Q E C» ((C Q (T) n U ) x (0,oo),^ >e (T)) , 

where Fj >e :— h J ^~ a (T) denotes the space of even, zero-mean functions with /i 2j+Q - 
regularity over T, j = 0, 1. We observe that Q(-, A) preserves symmetry about the 
y-axis, a property that follows quickly from the definition (2.3) of the governing 
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operator G and the representation of -0* derived in Remarks |5.1[ f). The fact that 
G preserves even functions is easily noted by the fact that every term of G(r) is 
either constant on T due to integration, or depends only on r, r XXl and r x , which 
are each even functions on T when r is itself taken to be even on T. 

Theorem 6.1. Fix £ G N. Then (0,£) is a bifurcation point for the equation 

(6.1) Q{f, A) = 0, (f, A) £ h 2 + e a (T) x (0, oo). 

In particular, there exists a positive constant 8i > and a nontrivial analytic curve 

(6.2) {(n(s),X e (s)) G h 2 + e a (T) xR: s G (-fc.fc), fa(0),A,(0)) = (0,*)} , 
such that 

Q(fe(s), \i(s)) = /or aZZ se(— 



and aZ/ solutions of (6.1) m a neighborhood of (0,£) are either a trivial solution 



(0, A) or an element of the nontrivial curve (6.2). Moreover, if X € R + \ N, then 



(0, A) is not a bifurcation point for (6.1 ). We can further conclude that 

a) lifting the curve (6.2), viaip+, we get elements of the family of 2tt/£ -periodic 
undulary curves with parameter values \B\ < 6i for some Se > 0. 

b) the bifurcation is a subcritical pitchfork type bifurcation. More precisely, for 
all £ G N, we have 

Xt(0) = and Xi(0) < 0, 
where " ' " denotes differentiation with respect to the parameter s. 



c) the bifurcating branch of undulary curves are unstable equilibria of (2.4), at 
least for parameter values \B\ < 8g sufficiently small. 

Proof. By working in the setting of even functions on T, we can take advantage of 
the cosine-Fourier series representation 

r(x) = f(k) cos(kz), for all f G /lo.eOOi 
fc>i 

where f(k) — — J T f(x) cos(kx)dx are the cosine-Fourier coefficients of f G Hq e (T), 



k > 1. It follows easily from (5.12) that the linearization DiQ(0, A) is a Fourier 
multiplier with symbol 

(M fc ) fe > 1 :=(A 2 -A; 2 )^ 1 . 
We then see that DiQ(0, A) is bijective from Fi j6 onto F 0y£ when A G M+ \ N, so 



that only points of the form (0,£) can possibly be bifurcation points for (6.1 ). 

Proceeding to verify that (0, £) is indeed a bifurcation point, we compute the 
kernel and range of DiQ(0,£) as 

Ni = span{cos(£c)} and Ri = span{cos(/cx) : k > 1, k ^ £}, 

respectively. Since h a (T) L2(T), we can borrow the L2-inner product to realize 
Ni as a topological complement to Ri as subspaces of Fi, e . By compactness of the 
resolvent R(X) := (A - £)G(r Jr ))~ 1 , A G p(DG(r*)), it follows that D X Q{Q,£) is a 
Fredholm operator of index zero. Then, defining the element vq :— cos(£x), and 
noting that 

Dl 2 g(0,£)v = -2£P v = -2^cos(fe) g i?^, 



it follows by [H Theorem 1.7], or HH Theorem 1.5.1], that (JeTT]) bifurcates at (0,£), 
which proves the main part of the theorem. 
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(a) Notice that ^(0, £) is precisely the constant function r* = l/£, and, utilizing 
the explicit characterization of equilibria from Section[4]above, we know that a curve 



of non-trivial equilibria for (2.4 1 containing this function must coincide with the 
family of ^-periodic unduloids. Hence, it suffices to show that Q(f, A) = if and 
only if G('0*(^, A -1 )) = 0, i.e. tp*(r, A -1 ) is an equilibrium of ( |2.4[ ) if and only if 
(r, A) solves (|6.1[). However, this follows immediately from the relation 



(6.3) G(r) = DMP r,r,)P G(r), r £ Ex, 



which was justified in deriving (5.11 ) above. 

(b) We will follow the characterization of bifurcation types as developed in [2T1 
Sections 1.6 and 1.7]. Computing the second derivative 

D\ s Q{Q t £)[% i %] = P (2£ 4 Q (cos 2 (£x)) - 2£ 3 P (cos 2 (£x)) - ^ 4 Q (sin 2 (fe)) 



- £ 3 sin 2 (£x) + 2£ 3 Q (cos 2 {£x)) 

e Po ( c ^m±l ) = _ £ 3 /2 cos{2ex)j 



we note that —£ 3 /2 cos(2£x) £ Re, from which it follows that A(0) = 0. Meanwhile, 
utilizing the representation [3TJ (1.6.11) and (1.6.9)] and following a considerable 
amount of computation, one will see that 

A(0) = -i^M = -/»<0, 
3 {2£vo,v } 

where (•, •) denotes the inner product on L 2 (T), within which Kq e (T) is embedded. 

(c) We can now track the so-called critical eigenvalue /x^(A) of the lineariza- 
tion DiQ{§, A), which is the eigenvalue which passes through at A = £, with 
non-vanishing speed, as guaranteed by the observed bifurcation. Moreover, via 
eigenvalue perturbation techniques, we can also track the perturbed eigenvalues 
/t^(s) of the linearization DiQ(fg(s), Xe (s)). Then, by taking a derivative of the 
relation (16. 31), we observe that 



(6.4) DG(Mr, v))DiMr, rj)h = Di^*(f , v)DiG(f, rj)h, h £ Fx, 

from which one will easily conclude that p,(s) must also be an eigenvalue of the 
linearization DG{^p^{fi (s), \J l {s)). Finally, by the subcritical structure of the bi- 
furcation, we conclude that, for sufficiently small |s| < 8e, the perturbed eigenvalue 
fie(s) has positive real part, from which instability follows by a similar argument 
to Theorem 15.41 above. □ 
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